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I.   INTRODUCTION 

In  any  inventory  model,  one  of  the  main  objectives  is 
to  avoid  both  overstock  and  out-of-stock  situations.   The 
latter  is  especially  critical  in  military  supply.   In  lot 
size  -  reorder  point  models  with  stochastic  demands,  a 
decision  has  to  be  made  as  to  when  to  reorder  and  how  much 
the  reorder  quantity  should  be.   The  goal,  then,  is  to 
achieve  in  some  sense  an  optimal  reorder  level  r  and  an 
optimal  reorder  quantity  Q.   One  way  to  accomplish  this  goal 
is  to  define  a  cost  function  K   which  then  acts  as  an  objective 
function  to  be  minimized  by  choice  of  Q  and  r  as  variables. 

Hadley  and  Whitin  [Ref .  1]  present  several  lot  size  - 
reorder  point  models.   The  one  of  interest  in  this  paper  is 
called  a  backorder  model  wherein,  if  demands  occur  when  the 
system  is  out  of  stock,  they  will  be  backordered  with  some 
penalty  cost,  to  be  called  the  backorder  cost .   Customers 
then  must  wait  until  all  orders  are  delivered  to  the  inventory 
manager.   Associated  with  delivery  is  a  delay  in  time  called 
lead  time .   In  this  paper  we  consider  only  the  case  where 
lead-time  demand  is  assumed  to  be  normally  distributed. 

Mathematically,  the  backorder  model  of  interest  here  has 
the  objective  function, 

K(Q,r)  =  ~  +   IC[§  +  r-vj  +  ^[(w-r)*(^-)  +  a<J>(^)]     (D 


where 

A  is  the  set-up  cost  (cost  of  placing  an  order) 

it  is  the  backorder  cost  per  unit  backordered 

X  is  the  average  annual  demand 

I  is  the  inventory  carrying  charge 

C  is  the  unit  cost  of  the  item 

u  and  a  are  the  mean  and  standard  deviation  of  the 

normally  distributed  lead-time  demand  random 
variable 

<t>(z)  is  the  standard  normal  density  function 

$(z)  is  its  complementary  cumulative  distribution  function, 

oo 

i.e.,   <Hz)  =  /   <J)  (x)  dx  . 
z 

K(Q,r)    is  called  the  average  annual  cost  and  is  composed 
of  three  different  types  of  cost:   the  ordering  cost,  the 
cost  of  carrying  inventory  and  the  cost  of  backorders.   These 
costs  are  represented,  respectively,  by  the  first,  second  and 
third  term  in  the  expression  for  K(Q,r)  in  (1). 

For  ease  in  the  algebra,  (1)  will  be  rewritten  as 


MQ,z)  =  ^  +  IC[§  +  az]  +  2A£  [(J)(2)  -  z»(z)]  (2) 


v.         r-u 
where   z  =  — -  . 
a 

Hadley  and  Whitin  [Ref.  1,  pp.  165  ff.]  point  out  that 

the  model  just  defined  is  appropriate  when  the  expected  number 

of  backorders  is  negligible.   They  then  propose  a  solution  to 

the  optimization  problem  by  claiming  that  K  is  convex  in  Q 


and  r  (or  z)  and  therefore  any  solution  obtained  by  setting 
the  partial  derivatives  of  K   equal  to  zero  will  determine  a 
global  minimum.   They  provide  an  iterative  scheme  for  such 
a  solution  along  with  the  assertion  that  the  solution,  when 
it  exists,  is  unique.   Finally,  they  discuss  the  fact  that 
a  solution  will  not  always  exist,  pointing  out  that  this 
"anomaly"  arises  because  a  backorders  term  was  omitted  in  the 
expression  for  the  carrying  cost.   Since  the  model  is  to  apply 
only  when  such  a  term  is  negligible,  they  assert  that  this 
situation  will  never  occur  when  high  backorder  costs  are 
involved  without  being  explicit  as  to  what  "high"  means. 

To  verify  these  results,  the  reader  is  directed  to  solve 
a  series  of  exercises,  the  first  of  which  is  to  show  that 
the  function  involving  the  backorder  cost  (the  third  term) 
in  (1) ,  say  J(Q,r) ,  is  convex  in  Q  and  r.   It  would  then 
follow  immediately  that  K  is  convex  in  Q  and  r.   But  it  was 
pointed  out  as  early  as  1964  by  Veinott  [Ref.  2]  that  J  is 
not  convex.   In  1969,  Brooks  and  Lu  [Ref.  3]  addressed  the 
same  problem  and  derived  a  general  result  which,  when  applied 
to  the  normal  case,  shows  that  J  is  convex  for  all   r  _>  y  , 
the  mean  lead-time  demand. 

However,  these  results  still  leave  the  main  issue 
unresolved.   In  the  first  place,  solutions  can  be  obtained 
where   r  <  y.    More  importantly,  the  real  issue  is  whether 


The  difference   r-y   is  sometimes  referred  to  as  a 
safety  level.   In  that  case,  a  negative  safety  level  simply 
means  it  is  more  economical  to  tolerate  a  backorder  position 
on  the  average  than  to  protect  against  that  position. 


or  not  K  is  convex,  and,  in  any  case,  how  its  minimum  value 
is  characterized.  The  purpose  of  this  thesis  is  to  clarify 
this  issue  and  characterize  the  solutions  to  the  problem  in 
a  more  explicit  manner. 

Section  II  will  show  that  in  fact  K    is  not  convex  and, 
in  addition,  will  exhibit  an  example  in  which  the  first 
partial  derivatives  -^r     and  ■*—  vanish  at  two  different 

oQ  dZ 

points.   Thus,  the  solution  is  not  unique  as  claimed  by 
Hadley  and  Whitin.   Section  III  will  deal  with  the  character- 
ization of  the  solutions  and  will  show  that  under  suitable 
conditions  there  are  always  two  distinct  solutions,  one  of 
them  being  a  minimum  and  the  other  a  saddle  point.   Conditions 
for  the  existence  of  a  solution  will  be  given  in  terms  of  the 
set-up  cost  A.   Section  IV  will  discuss  the  application  of 
the  results  in  Section  III  to  practical  inventory  problems 
along  with  the  iterative  scheme  given  by  Hadley  and  Whitin  to 
solve  for  a  solution. 


II.   COUNTEREXAMPLES 

As  pointed  out  in  the  introduction,  the  backorders  term 
in  the  cost  function  is  not  convex  as  asserted  by  Hadley 
and  Whitin  [Ref.  1,  Problem  4-6,  p.  221],   Even  so,  it  might 
still  be  the  case  that  K   is  convex  as  asserted  by  the  authors. 
That  this  is  not  so  will  be  demonstrated  in  part  A.   On  the 
other  hand,  even  though  K   is  not  convex,  it  might  still  be 
the  case  that  the  partial  derivatives  vanish  at  a  unique 
point  and  that  point  determines  a  minimum.   That  this  is  not 
so  will  be  shown  in  part  B  and  a  characterization  of  solutions 
will  be  given  in  the  next  section. 


A.   CONVEXITY  CONDITION  FOR  K 

The  function  K.    is  convex  if  and  only  if  its  Hessian 
matrix  Q   is  positive  semi-definite  [Ref.  4],   Since  K  is  a 
function  of  two  variables,  the  matrix  Q   will  be  positive 
semi-definite  when  its  diagonal  terms  and  its  determinant 
are  nonnegative. 

For  the  case  at  hand,  using  (2) ,  the  Hessian  matrix  for  K 
is  given  by 


n  = 


~    (A  +    Trn(z)) 


irA<J>  (z) 


ttA$  (z) 


Q' 


TT  X  4>  (z) 

cQ 
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where    n"(z)  =  a<|>(z)  -  az$(z) 

or       "n(z)  =  au(z)  (3) 

where    u(z)=<J>(z)-z$(z)   .  (4) 

In  Appendix  A  it  is  shown  that   u(z)  >  0   for  all  z  . 
Hence,  for  nonnegative  set-up  costs  A,  it  is  clear  that  the 
diagonal  terms  of  ft  are  positive.   Thus  K    is  convex  if  and 
only  if   det  Q     _>   0  . 

Now, 


det  fl  =  21TX2^Z)  (A  +  TrrT(z))  -  »*2»*<«> 

aQ  Q4 


and  thus   det  R  >  0   if  and  only  if 


2**VZ>  (a  +  Tin(z))  -  ^2xV(z)  >  o 

aQ4  Q4 

4  ... 

Multiplying  the  above  inequality  by   Q   >  0   and  dividing 

2 

by      ttX      >    0      gives 


^^-)  (A  +    irn  (z)  )    -    tt<D2  (z)     >      0 


PA  9  9 

or  ~  <}>(z)    +   2tt<J>    (z)    -    2ttzcj)(zH(z)    -   tt<T(z)    >      0  (5) 


Defining  V(z)  to  be  equal  to  the  left-hand  side  of  (5),  then 
the  following  summarizes  convexity. 


11 


RESULT:   The  function  K    is  convex  if  and  only  if   det  Q    >_   0 
or  equivalently  if  and  only  if  V(z)  ^L  0  • 

REMARK:   This  result  holds  for  any  Q  >  0  and  it  is  shown 
in  Appendix  A  that  V(z)  >  0  whenever  z  >_  0  . 

A  counterexample  to  the  convexity  of  K   is  then  attainable 
from  an  example  provided  by  Hadley  and  Whitin  [Ref.  1,  p.  173]. 
There  the  parameters  chosen  are  t\   =  2000,   A  =  1600  ,   C  =  50, 
I  =  .20  ,   A  =  4000  ,   y  =  750  ,  a   =  50  .   For  such  a  case 
V(z)  becomes 


V(z)  =  1604>(z)  +  4000cf>2(z)  -  4000z<J>(z)$(z)  ~  2000$2(z) 


and  for   z  =  -1   (equivalently,  r  =  700) 

V(-l)  =  38.72  +  234.256  +  814.3784  -  1415.57138 
=  -328.21698  <  0 

Thus,  K    cannot  be  convex. 

B.   COUNTEREXAMPLE  TO  UNIQUENESS  OF  THE  SOLUTION 
DEFINITION.   (Q,z)  will  be  called  a  solution  whenever  Q  >  0 
and  z  satisfy  the  following  equations: 


„    8K      XA    IC    7TXn(z)      AA  IC    TTAau(z)          ,cx 

0  ■  3q  =  -  ?  +  -  -  —j-  •  -  ?  -  "  -—            <6> 

0  -  H=  (IC  -  »*»<»> )0  (7) 

dZ  Q 
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It   is    convenient   to  write    (6)    and    (7)    as 


t|2X(A    +    TTTjJz))        _      tj2A  (A    +    ttgu(z)  )  ,0. 

Q    =     y  IC  "      y  IC  (8) 

and 

•(>)    =  ^  (9) 

COUNTEREXAMPLE.   Let  tt  =  A  =  a  =  1  ,   I  =  .1  ,   C  =  2.42 

1        2 

and  y  =  3.   Let  A  =  2 (  1) (2  42)  $  (-05)  -  u(.05)  =  .1018  , 

where  u(z)  is  defined  in  (4) .   Then  there  are  two  distinct 

solutions 

•<z.) 
1)   (Q1,z1)  =  (  (>l)(2t42),  zx)   where   Z]_  =  .05 


*(z  2) 
2)   (Q2'z2)  =  (  (  x) (2  42) '  Z2*   where   -1.81  <  z2  <  -1.80 


To  see  that  (0-,z,)  is  a  solution, 


2 

1   r    $   (.05)  ,  n(-N  -i 

L2(.l)(2.42)   '  U^-UD'J  +  (.1)(2.42)    (1)  (1)  (l)u(.05) 


Q1,z1  _$2(.05)  $2(.05) 


(.1)2(2.42)2  (.1)2(2.42)2 

(.1)(2.42)  .  (.l)2(2.42)2u(.05)    (.1)(2.42)    (.l)2  (2.42)  2u(  .05) 
_ +  _ 1 + 


$2(.05)  2  S>2(.05) 


=  0  .   Thus,   equation  (6)  is  satisfied. 
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3K 
Tz" 


=  (.1X2.42)  -  (1)ffio5)°5)  =  (.D(2.42)  -  (.1)  (2.42)  =0 
Q1'Z1  (.1)  (2.42) 


Hence  equation  (7)  is  satisfied.   Therefore,  (Q,  ,z,)  is  a 

solution. 

1       2 
Now,  let   S(z)  =  2 (  l) (2    42)  $  ^  ~  u^z^  ■   Since  S  is 

clearly  a  continuous  function  and 


S(-1.81)  <   .0997  <   A  <  .1060  <  S(-1.80) 


there  must  exist  some  z2  such  that   -1.81  <   z2  <  -1.80 
and  S(z2)  -  A.   The  verification  for  (Q2,z2)  is  done  in  a 
similar  manner. 


M 

9Q 


$   (z2) 
1    [2(.l)  (2.42)     -u(z2)]     +    (.1)(2.42)         (1)  (1)  (l)u(z2) 


Q2/Z2  _     $2(z2)  2  $2(z2) 


(.1)2(2.42)2  (.1)2(2.42)2 


(.1)  (2.42)   +    (.l)2(2.42)2u(z2)   +    (.1)  (2.42)        (.I)2  (2.42)2u(z2) 


$2(z2)  2  $2(z2) 


9K 
8z 


=  0    .      Thus,    equation    (6)    is    satisfied. 


=   (.1)    2.42     -   (  }Bz?\  =     .1     2.42     -      .1     2.42     =0 

O     7  ^   2) 

U2'z2  (.1)(2.42) 


and  hence  equation  (7)  is  satisfied.   Therefore,  (Qo'z2^  ^S 
also  a  solution . 

14 


To  illustrate  these  results,  a  graph  similar  to  the  one 
given  on  p.  171  of  Ref.  1  is  reproduced  in  Fig.  2-1  for  this 
special  case.   In  this  graph   r  =  az  +  y   is  expressed  as  a 
function  of  Q  for  each  of  Equations  (8)  and  (9) .   It  shows 
that  the  two  curves  intersect  at  more  than  one  point. 


Fig.  2-1. 
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1        2 

The  value   A  =  of  lW2  A2)     ®    ^«05)  ~  u(.05)   was  by  no 

means  chosen  by  accident.   As  it  will  be  seen  in  the  next 
section,  the  set-up  cost  A  plays  an  important  role  in  the 
characterization  of  the  solutions,  and  it  will  be  shown  that, 
under  appropriate  conditions  of  the  other  parameters, there 
are  always  two  distinct  solutions,  one  of  which  is  a  minimum. 
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III.   CHARACTERIZATION  OF  MINIMUM  SOLUTION 

Let  a  set  of  parameters  be  given.  To  simplify  matters, 
define 


S(z)  =  —^   <D2(z)  -  ttou(z)  (10) 


A  series  of  lemmas  are  convenient  to  derive  the  main  result 
in  a  theorem. 

LEMMA  1.   Let  (Q*,z*)  be  a  solution.   Then, 

(a)  A  -  S(z*) 

(b)  V(z*)  >  0  if  and  only  if  *(z*)  >  ^S-   and  V(z*)  =  0 

if  and  only  if  <t>(z*)  =  — r-  . 

T  TT  A 

Proof  (a).   By  definition  of  a  solution,  Q*  and  z*  must 
satisfy  (8)  and  (9)  of  Section  III.   Thus, 


2TT(A+7Tau(z*))    fn4,2    (ttA)  $2(z*) 

(Q*) 

1C  (IC)^ 


or 


2 

2AA  +  27rXau(z*)  =  (ir^   $2(z*) 


2AA  =  (7T^   $2(z*)  -  27iAau(z*) 


2 
A  =  ~£  *  (z*)  -  7rau(z*)  =  S(z*) 
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Proof  (b)  . 


?A  ? 

V(Z*)  =  ~    <$>(Z*)     +  27T({>(Z*)  (*(2*)-Z**(z*))  -  7T$   (Z*) 


?A  2 

■^  (})(Z*)   +  27T(J)(z*)u(z*)   -  7T<r(Z*) 


=  2(J»(Z*){|+  7TU(Z*))   -  TT$2(Z*) 

2 
=  2((>(z*)  (J^$2(z*)-ttu(z*)+ttu(z*)  )  -  7T$2(z*)  from  part  (a) 

=  |~  (^(Z*)$2(Z*)   -  TT$2(Z*) 

=  ^*2(z*)  [<Kz*)  -  *£] 

7T  /    2 

Since  Yfv7  $  (z*)  >  ^,  the  result  stated  in  (b)  follows 
immediately. 

Q.E.D. 
It  is  not  always  true  that  there  exists  a  solution 
for  any  given  set  of  parameters.   The  following  lemma  discusses 
conditions  for  which  K.   cannot  be  minimized. 

LEMMA  2.   Given  A  >_  0  and  it,    A,  I,  C,  o  all  positive,  a 
necessary  condition  for  the  function  K   to  be  minimized  is 
that  2££  <  (f)(0)  . 

Proof.   The  proof  is  divided  into  two  cases. 

ICo 
Case  (i)  :   Suppose  — r-  >  4>(0)  and  K  is  minimized  at  (Q*,z*). 

Then  the  partial  derivatives  must  vanish  at  (Q*,z*).   But 
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2. 
V(z*)  =  — ^  *2(z*)  [<Kz*)  ~  =^]  <  0 


since  (f)(z*)  _<  <J>(0)  and  so  (Q*,z*)  is  a  saddle  point.   This 

is  a  contradiction. 

ICcr 
Case  (ii)  :   Suppose  — r-  =  <£>(0)  and  K   is  minimized  at  (Q*,z*) 

Again,  the  partial  derivatives  must  vanish  at  (Q*,z*).   Now 

if  z*  j-   0  then  V(z*)  <  0  and  hence  (Q*,z*)  is  a  saddle  point 

Thus  it  must  be  the  case  that  z*  =  0.   But  then,  by  Lemma  1, 

part  (a) , 


A  =  hz  <J,2(0)  "  Trau(0) 


2, 


=  Sic  "  ^^ 


TTQ       TTX  .... 

.     -  „„♦«» 


TO  1  . / A\ 


8    <j>(0) 


,1  -  8d)2(0) 

7Ta( ^Tol ) 


2 

Since  1-84  (0)  <  0,  A  <  0  contradicting  the  fact  that 


A  >  0. 


ICa 
Thus  K   cannot  be  minimized  unless  — y  <  <f>(0) 


Q.E.D. 
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LEMMA    3.      Suppose      ~^-   <    (J>(0).      Let      z^    >    0      and      -z^    <    0 

ICO 

be  the  unique  points  such  that  $(-z    )    =    ^>(z)    =  — r-  . 

O  O       7T  A 

Then,  with  S(z)  defined  in  (10), 

(a)  S  has  a  minimum  at  z  where  A'  =  S(z  )  <  0 

(b)  S  has  a  maximum  at  -z  where   AQ  =  S(-z  )  >  0 

(c)  S  is  strictly  concave  increasing  on  the  interval  (-°°,-z  ) 


Proof.   S(z)  =  J^  $2(z)  -  7rau(z) 


7T2X 
S"  (Z)  =  -  ~^  <f>(z)$(z)   -  TTQU'  (Z) 


IC 


<J)(z)$(z)  -  TTcr  (-$  (z)  )    (see  Appendix  A) 


TT2X 
=  TTO$(z)   -  -j^-  (|>(z)$(z) 


2 

7T  X  .  ,  .  rICa    ,  .  .  . 
$(z)  [— t 4»(z)  ] 


IC     v  '   l  TT, 

1T2X 

Since  -j^r   §(z)  >  0  ,   the  first  derivative  S'(z)  vanishes 

only  at  -z   and  z  .  Also,  S'(z)  <  0   for   -z   <  z  <  z0   and 

S1 (z)  >  0   for   z  <  -z   or   z  >  z  ;   that  is,  S  is  decreasing 

in  the  interval  (-z  ,z  )  and  increasing  in  the  intervals 
(-°°,-zQ)  and  (zQf«>)  . 
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2  2 

S"(z)    =  -  ?L±  <|>(z)  fiSjL  -   <j>(z))    +  !L£  »(Z)  (z4»(z)) 


=  -ttq^z)    +   e!A|1^J    +   ^2Xz4)(z)Mz) 


2 

U     ^(Z)       (<|>(Z)      +     Z$(z))      -     TTOCHZ) 


ira<J>(z)  [^    (<{>(z)    +   z*(z))    -   1] 


Thus,    at   zQ, 


S"(z0)    =   ^*(z0)[^    (JC£+    z0*(Zq))    -   1] 


«a#(«0)[l  +  jgL  z0*(z0)    -   1] 


T^  zo*(zo)0(zo)    >      0  (11) 


Therefore.  S  has  a  minimum  at  z  . 

o 

To  show  that  A^  =  S(zQ)  <  0  ,   suppose,  to  the  contrary, 
that   S(zQ)  >  0  .   Then,  since   S'(z)  >  0   for  z  >  z  , 
S(z)  >  S(zQ)  >  0   for  all   z  >  zQ   and  hence   lim  S(z)  >  0. 

Z->oo 

But   lim  S(z)  =   lim  ~F   <J>2(z)   -  -no   lim  u(z)   =   0   (see 
Z-*"00  Z-*"00  z-*00 

Appendix  A).   Thus,   A'  =  S(z  )  <  0  .   This  together  with 

(11)  proves  part  (a) . 
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Now  at  -zn, 
o 


S"(-z0)  =  w#{-«0)[^  (££-  z0*(-z0))  -  1] 


=  w#(-«0)[l  -  ^  •0»(-«o»  -  x' 


z„*(-z„)*(-z„)  <  0  .  (12) 


IC  "oT  x  "o' " x  "o 
Therefore,  S  has  a  maximum  at  -z  . 

To  show  that  A  =  S(-z  )  >^  0  ,   again  suppose  that 

S(-zQ)  <  0.   Since   lim  S(z)  =0  ,   there  exists  some  z1 

z-*-oo 

such  that   S(z')  >  S  (-z  )  ,  for  otherwise   S(z)  <_   S  (-z  )   for 
all  z  implies  that   lim  S(z)  <_  S(-z  )  <  0  .   But  the  existence 
of  such  z1  contradicts  the  fact  that  S  is  maximized  at  -z  . 
Thus,   A   =  S(-z  )  >  0  .   This  result  coupled  with  (12)  proves 
part  (b) . 

Part  (c)  can  now  be  easily  proved  by  noting  that  for 
z  <  -zQ  <  0  , 


4>(z)  +  z$(z)  <  4>(z)  <  <()(-zo)  =  ~2 


hence  ~L_  ($(z)  +  z$(z))  <  1 


or     j~    (<j)(z)  +  z$(z))  -  1  <  0 


and   therefore      S"  (z)    =    tto6  (z)  [^-    (d>(z)    +    z$(z))    -    1]    <    0 

r     iCa 

Thus,  S  is  strictly  concave  increasing  on  the  interval  (-°°,-z  ) 

The  proof  of  Lemma  3  is  now  complete. 

Q.E.D. 
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REMARK ;   The  minimum  and  maximum  in  Lemma  2  are  clearly  unique. 
Moreover,   lim   S(z)  =  -00   (see  Appendix  A).   A  graph  of  S(z) 


z-*—( 


corresponding  to  the  example  of  Section  II  is  given  in  Fig.  3-1 


Fig.  3-1. 

The  following  lemma  shows  the  importance  of  A  in  the 
characterization  of  a  solution. 


LEMMA  4 .   Suppose 


ICa 


<  $(0)   and  A  >  0   is  given.   Then 

TTA  —  3 

(Q*,z*)  is  a  solution  if  and  only  if  A  =  S(z*). 

Proof.   The  necessary  part  has  been  proved  by  Lemma  1,  part 

(a) .   The  sufficient  condition  will  now  be  proved.   Suppose 

2 
A  =  S(z*)  -  j-£   $2(z*)  -  ttou(z*)  .   Let   Q*  =  l*^z*>  . 

Then,  it  can  be  easily  shown  by  substitution  of  A  and  Q* 
into  equations  (6)  and  (7)  that  (Q*,z*)  is  a  solution. 

Q.E.D. 
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REMARK;   It  should  be  observed  immediately  from  Lemmas  4  and  3 

that  a  necessary  condition  for  a  solution  to  exist  is  that 

0  _<  A  £  A  .   The  following  theorem  summarizes  all  the  results. 

THEOREM.   Given  tt,  X,    I,  C,  a  all  positive  and  — £  <  <f>(0) 

7T  A      T 

and  A   >  0  ,   then  for   0  <  A  <  A   there  exist  exactly  two 
solutions  (Q-,*,z,*)  and  (Q2*,z2*)  such  that   z,*  >  z2*  ; 
(Q-,*/Z,*)  is  a  minimum  and   (Q2*,z2*)   is  a  saddle  point  for  K. 

Proof.   Since  S(z)  is  strictly  concave  increasing  in  the 
interval  (-oo/-z  )  ,  and  strictly  decreasing  in  the  interval 
(-z  ,z  )   with   lim  S(z)  =  -°°  ,   S  (-zQ)  =  AQ  >  0   and 
S(z  )  =  A'  <  0  ,   the  function  S(z)  must  have  two  roots  z-, 
and  z2  where   z,  <  -z    and   -z   <  z2  <  z   .   Thus,  for  any 
A  such  that   0  <_  A  <  A   ,   there  are  two  distinct  values 
z-j*  and  z2*   such  that   z-,  <  z2*  <  -z   <   z-^*  <  z2   and 


A  =  S(z  *)  =  S(z2*).   By  merely  defining  Q1*  = 


ttA$  (z,  *) 
IC 


ttX<J)  (z2*) 
and  Q2*  =  « '      ■*•*  follows  by  Lemma  4  that  (Q,*,z,*) 

and  (Q*,z2*)  are  solutions.   But,   -z   <  z-,*  <  z    implies 

that   <j>(z  *)  >  ^p,   and  hence  by  Lemma  1,  part  (b)  ,  V(z1*)  >  0 

Thus,  (Q-j*,z-|*)  minimizes  K .   Similarly,   z2*  <  -z    implies 

that   d>(z~*)  <  — ?-  and  hence  by  the  same  lemma  V(z0*)  <  0  . 

Thus,  (Q9*/Z*)  is  a  saddle  point. 

*  Q.E.D. 

REMARK:   For   A  =  A   =  S  (-z  )  ,  the  only  solution  possible  is 
determined  by   z*  =  -zQ  .   But  then,   <}>(z*)  =  -—■  and  thus 
V(z*)  =  0 ,  so  that  no  conclusion  can  be  drawn. 
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IV.   COMPUTATIONAL  ASPECTS 

Before  attempting  to  solve  a  given  problem  involving  a 
set  of  parameters,  some  simple  tests  should  be  performed  to 
see  whether  or  not  there  is  any  solution  ~ 

1)  The  first  test  is  ^£  <  d>(0). 

2)  If  this  test  is  passed,  then  find  -z   and  A  =  S(-z  ). 

Two  cases  could  then  occur. 
Case  (i) ;   AQ  =  0.   The  only  solution  then  is  -z   and  V(-z_)  =0 

and  it  has  been  pointed  out  that  no  conclusion  may  then  be 
drawn . 

Case  (ii)  :   A   >  0.   If   0  <_  A  <  A  ,  then  there  exist  exactly 
two  solutions  (Q,*,z1*)  and  (Q2*/Z2*)  with   z2*  <   "zo  <  zl*  * 
As  shown,   (0*,z,*)  is  a  minimum  and  (Q2*/z2*)  ^"s  a  Ea(^dle 
point.   If  A  =  A  ,  then  again  the  only  solution  is  -z   and 
no  conclusion  can  be  drawn  as  in  case  (i) . 

Hadley  and  Whitin  give  an  iterative  scheme  for  solving 
for  (Q*,r*)  in  Ref.  1.   This  scheme  seems  to  converge  always 
to  the  minimum  value  and  never  to  the  saddle  point.   But  there 
is  a  simple  explanation  for  this.  The  scheme  begins  by  choosing 

initially  a  Wilson  0_w  -  \J^Yc~     ^or  t^le  va-^ue  °^  °-  to  ^e 
used  in  Eq.  (9)  (see  Fig.  2-1)  and  Q„  <  O* .   The  iteration 
then  proceeds  to  choose  a  value  of  r  corresponding  to  Q„  from 
Eq.  (8) .   With  that  value  of  r  a  new  value  of  Q  is  selected 
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from  Eq.  (9)  and  it  will  always  be  the  case  that   Q  <  Q,* 
in  that  selection.   Proceeding  in  this  manner,  the  pairs 
(Q,r)  converge  to  the  solution  (Q-,*,r.,*)  and  cannot  yield 
(Q2*/^2*^*   Thus,  when  a  minimum  exists,  this  computational 
scheme  may  be  relied  upon. 

However,  in  any  case  of  doubt,  a  simple  test  is  furnished 

ICg 
by  checking  whether  or  not   cf>(z*)  >  -=j~     ^or  t*ie  solution. 

It  has  been  pointed  out  that  the  boundary  case  A  =  A 

causes  problems  in  that  no  conclusion  may  then  be  drawn.   This 

presents  no  difficulty  for  practical  uses,  however.   After 

all,  an  optimal  solution  may  be  given  for  any  A  arbitrarily 

close  to  A  .   In  an  application  it  would  be  a  rare  circumstance 
o 

that  produced  a  set  of  parameters  in  which  A  (typically 
measured  in  dollars  and  cents)  would  be  exactly  A  . 
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V.   CONCLUSION 

While  the  results  of  this  thesis  only  apply  to  the  case 
of  normal  lead-time  demand,  they  are  of  more  than  mere 
academic  interest.   In  fact,  that  is  a  common  assumption  made 
by  Hadley  and  Whitin  along  with  users  of  the  formulas  and 
schemes  given  by  them.   And  it  is  a  happy  coincidence  that 
the  algorithms  given  are  indeed  valid  even  though  the 
assumptions  under  which  they  were  derived  were  replete  with 
errors . 

Nevertheless,  it  is  a  recommendations  of  this  thesis  that 
further  study  into  other  types  of  lead-time  distributions  be 
made.   One  of  the  key  features  of  the  results  derived  here 
was  the  examination  of  the  Hessian  matrix  with  the  observation 
that  its  positivity  did  not  depend  on  Q.   This  will  still  be 
true  with  a  density  other  than  normal  and  so  much  of  the  same 
general  technique  ought  to  apply. 

One  of  the  features  of  this  thesis  has  been  to  amplify 
and  clarify  completely  the  minimum  solution.   Conditions  are 
given  for  existence, and  formulas  are  given  along  with  simple  test 
criteria  when  existence  is  assured.   A  by-product  of  the 
results  is  the  fact  that  negative  safety  levels  do  indeed 
occur  as  the  unique  optimal  solution.   In  fact,  a  simple 
examination  of  a  graph  like  S(z)  will  show  immediately  and 
simplify  those  values  of  the  set-up  cost  A  for  which  such 
solutions  are  valid.   Since  the  solution  in  such  cases  is 
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unique,  certainly  negative  safety  levels  should  not  be  ignored 
or  tampered  with  (truncation  to  \i    for  example)  in  an  application 
It  should  also  be  noted  that  the  solutions  are  characterized 
independently  of  u. 

As  a  final  recommendation,  this  thesis  was  restricted  to 
the  so-called  backorders  model  where  a  customer  must  wait  for 
delivery  when  the  inventory  position  is  out  of  stock.   There 
are,  however,  other  models  such  as  the  "lost  sale  case" 
discussed  by  Hadley  and  Whitin  to  apply  in  this  situation. 
An  analysis  similar  to  that  applied  here  is  surely  called  for 
so  that  resultant  solutions  are  validated  and  characterized. 
For  the  authors  continually,  appeal  to  convexity  for  drawing 
conclusions  in  those  cases. 
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APPENDIX    A.       VARIATION    OF    V(z) 


As    in   Section   II, 


P. 


~     (A    +     TTTI(Z)) 

Q 


pX$(z) 


ttA<I>(z) 
Q2 


ttA<j)  (z) 
OQ 


The   diagonal   term     — ^Qz        is   certainly  positive. 
rj"(z)    =   a(<f>(z)    -    z$(z))    =   cju(z)      as   defined   in    (4) 


u1  (z)    =  -z<f>(z)    -    (*(z)    -   z<J>(z)) 


=  -$(z)    <    0  for  all    z 


u"  (z)    =    4>(z)    >    0  for   all   z 


Clearly,      lim     u(z)    =    lim      (<{>(z)    -    z$(z))    =   0   +  °°  =   +^ 

2~>-— co  Z-*-00 


and 


lira     u(z) 

z-*-co 


lim    (4>(z) 

z->co 


-    z$(z) )    = 


-    lim  z$(z) 

2-vco 


-    lira 

z-*oo 


1 
$TzT 
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Applying  L1 Hospital's  Rule  repeatedly, 


lim  u (z)    = 

£->-0O 


-lim 
z-*-°° 


<f>(z) 

$2(z) 


=   -   lim 

z->oo 


$2(z) 
«J>(z) 


-    lim   -2Mz)f(Z) 
„_^«        -z4>(z) 


=  _  lim  mil  =  o  . 

z 

2-K» 


Thus,  u(z)  is  a  strictly  convex  decreasing  function 
approaching  zero  as   z  ■*■  +°°   and  hence   u(z)  >  0   for  all  z 
A  sketch  of  u(z)  is  given  in  Fig.  A-l. 


Fig.  A-l 


Since   u(z)  >  0   for  all  z,   n(z)  =  au(z)  >  0   for  all  z 

2  \  

and  therefore  the  diagonal  term  — «•  (A  +  ttti(z))  >  0   for 

Q 
all  z  since   A  >  0  . 
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2A  2  9 

V(Z)      =    ~    4)(Z)      +     2tT<J>      (z)      -     27TZ(f)(z)  $(Z)      -     TT$^(z) 


Define      T(z)    =    2(})2(z)    -    2z<t>(z)<Mz)    -    $2(z)     ,      then 


?A 
V(z)    =   ~   <J>(z)    +    ttT(z) 


T'  (z)    =    -4z(})2(z)    +    2z2({>(z)  $(z)    +    2z<f>2(z) 


=    2z2cf)(z)  <Mz)    -    2z<J>2(z) 


=    2z4>(z)  (z$(z)    -    4>(z)  ) 

=    -2z4>(z)u(z)  (A.l) 


T"  (z)    =   -2<J>(z)u(z)    +    2z2^(z)u(z)    +   2z<J>(z)$(z)  (A. 2) 


Thus, 


9A 

V  (z)    =   -  ~   z<j>(z)    +    ttT'  (z) 

2A 
=   -   —  z$(z)    -    2ttz<J)  (z)  u  (z)       from    (A.l) 

=   -2z<J>(z)  (    |  +    ttu(z)) 


Since      u(z)    >    0  for   all  z 

>    0  for  z  <  0 

=    0  for  z  =  0 

<    0  for  z  >  0 


31 


V"(z)  -  - 


2A         2A   2 

—  <f>(z)  +  —  z  <J>(z)  +  ttT"  (z) 


At   z  =  0  , 


?A 
V"  (0)  =  -  ~  <J,(0)  +  ttT"  (0) 


2A 


'—   <J>(0)  -  2tt(|>(0)u(0)     from  (A. 2) 


Thus   V" (0)  <  0   and  hence  V(z)  is  maximized  at  z = 0   and 
V(0)  =  •—  (f)(0)  +  ttT(O)  =  ^  (.3989)  +  tt(.0683)  >  0  . 


lim  V(z) 

Z-M-oo 


0  +  tt  lim  T(z) 

Z-H-oo 


=  0 


lim  V(z) 


0  +  tt  lim  T(z) 


=  TT(-l)   =  -  TT 


A  sketch  of  V(z)  is  given  in  Fig.  A-2 


Va 


Fig.  A-2. 

Since  V(z)  attains  a  maximum  with  positive  value  at  z  =  0  and 

is  decreasing  for  z  >  0  and   lim  V(z)  =0  ,   V(z)  >  0  for  all 

z->+°° 
z  >  0. 
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